Influence of the isospin and hypercharge chemical potentials on the location of the 
CEP in the fi B — T phase diagram of the SU(3) L x SU(3) R chiral quark model 
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We investigate the influence of the asymmetric quark matter (p u 7^ pd 7^ p s ) on the mass of the 
quasiparticles and the phase diagram of the chiral quark model parametrized at one-loop level of the 
renormalized theory, using the optimized perturbation theory for the resummation of the perturba- 
tive series. The effect of various chemical potentials introduced in the grand canonical ensemble is 
investigated with the method of relativistic many-body theory. The temperature dependence of the 
topological susceptibility is estimated with the help of the Witten-Veneziano mass formula. 
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INTRODUCTION 



The study of a system of particles at finite density and temperature is phenomenologically interesting because in 
heavy ion collision experiments the initial state is such that the chemical potentials /is, /!/, [iy (conjugate to the baryon 
charge, third component of the isospin and hypercharge, respectively) are non vanishing, although the last two are much 
smaller than the first one. Assuming thermal equilibrium, thermal models show that the strangeness chemical potential 
in central Si+Au collisions at the Brookhaven AGS experiment was 20-25% of the baryonic chemical potential for 
which the best fit gives fiB = 540 MeV [1] . For central Pb+Pb collisions at CERN SPS experiments the value of the 
strangeness chemical potential was ~ 25 — 30% and that of the isospin chemical potential ~ 2 — 5% of the value of 
Hb estimated to be around 233 — 266 MeV [2, 3]. The Compressed Baryonic Matter (CBM) experiment at FAIR in 
Darmstadt will explore regions of the QCD phase diagram with moderate temperature up to such high values of the 
^ ' baryonic density which are comparable with those in the core of neutron stars [4]. 

In many-body theory chemical potential is introduced to any conserved charge. In heavy ion collision experiments 
the baryon number, isospin and hypercharge can be considered conserved due to the short time elapsed between 
the formation of the fireball and its freeze-out, during which only the strong interactions play important role, the 
electroweak interactions being negligible. It is expected that in the very early stage of the fireball's evolution strangeness 
is abundantly produced in the deconfined phase through gluon-gluon fusion [5] , while in the hadronic phase in the vicinity 
of the transition multi-mesonic reactions will play an important role in the fast redistribution of strange quarks [6]. 

The influence of the isospin chemical potential on the chiral phase transition is currently actively investigated, because 
this effect can in principle be tested experimentally. As noticed in [7], using different isotopes of an element in heavy 
ion collision experiments will vary fij keeping fxs constant. Moreover, the system with real fij represents no extra 
difficulty in lattice field theory compared to the introduction of \xb- For two flavors the simulations at [ib =0 and 
fii 7^ is not even affected by the sign problem [8]. For (i; / a generic result coming from effective models of 
the strongly interacting matter without the U(1)a anomaly appeared to be the splitting in the [Ib — T plane of the 
first order transition line into two transition lines. This effect was observed in random matrix model [9], NJL model 
[10], strong coupling limit of the staggered lattice QCD [11], all with two flavors and in the three flavors ladder QCD 
[12]. This would imply the existence of not only the two phases having (in) 7^ 0, (dd) 7^ and (uu) = (dd) = 0, 
respectively, but also of a phase with (uu) = and (dd) 7^ 0. It was shown in [13, 14] that the structure with two 
transition lines and critical end points ceases to exist for a sufficiently strong U(1)a breaking, above which the two 
strongly coupled condensates vanish simultaneously. In a hadron resonance gas model it was found that at fixed baryon 
chemical potential the pseudocritical temperature of the transition between the hadronic and the quark-gluon plasma 
phases is lowered as either the isospin or the strangeness chemical potential is increased [15]. 

Because of their phenomenological implications, it is natural to study to what extent these results are present in 
another low energy effective model, the chiral quark model, widely used for studying the chiral behavior of strongly 
interacting matter. In the past few years we have investigated the thermodynamics of this model for two and three 
quark flavors at \xb = and /is 7^ 0, while leaving = (xy = [16-19]. As a continuation of these previous studies, 
in this paper we consider the influence of the chemical potentials on the chiral phase transition up to such high values 
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of the isospin chemical potential above which the condensation of pseudoscalar mesons occurs. The pion and kaon 
condensation phase, which is beyond the scope of our present investigation, was studied both with lattice methods and 
using effective theories [20-26]. 

The paper is organized as follows. In section II we present the model, its one-loop parametrization and the intro- 
duction of the chemical potentials. The variation of the location of the critical end point in presence of /i/ and /iy is 
studied in section III. There we investigate also the temperature and density dependence of the one-loop pole masses 
of the pseudoscalar mesons. We conclude in section IV. 



II. THE SU(3)l x SU{3)r SYMMETRIC CHIRAL QUARK MODEL 

The Lagrangian of the model containing explicit symmetry breaking terms is 

L = ^Tr (d^d^M + mlMHl) - f t (Tr (M^M)f - / 2 Tr (M f M) 2 - g (det(M) + det(M f )) 

+^(i^-g F M 5 )^. (1) 

The constituent quarks are contained in the field ip: tp = (u,d,s). The two 3x3 complex matrices are defined in 
terms of the scalar <7j and pseudoscalar 7Ti fields as M = ^= X)i=o( <T ' ; + an d ^5 = \ Ei=o( <Ji + il^i)Xi, with 

Xi : i = 1 ... 8 the Gell-Mann matrices and Ao := v/|l- The fields with well defined quantum numbers are obtained 
with a block-diagonal transformation f a = T ai fi, / g {a, 7r}, where T = diag(l, r, 1, r, t, 1) and r = f j f\. As 

a goes from to 8, the components of the scalar and pseudoscalar fields go trough ao, a>o, %, 03, k + , k~ , k° ,R° , a$ 
and 7To, 7T+, 7T _ , 7r 3 , K + , K~, K a , K°, ir s , respectively. The physical fields ir° (neutral pion), 77 and 77' mesons in the 
pseudoscalar sector and a[j (neutral a ), a and /o in the scalar sector are obtained as linear combinations of the 
corresponding fields in the two mixing 0,3,8 sectors. 

In this paper we investigate the pattern of symmetry breaking realized in nature, with the SU(3)a x U{1)a x SU (3)y 
symmetry completely broken, that is the isospin SU(2)y is also broken. In addition to the spontaneous symmetry 
breaking, explicit breaking is also considered with the introduction of external fields for all the diagonal generators of 
the scalar sector. This results in having three non- vanishing condensates in the broken symmetry phase: v$ = (as), for 
S = 0, 3, 8. The condensates determines the tree-level scalar and pseudoscalar masses: 

m S.ap = m2 &al3 ~ 6G a/ 3 7 U 7 + AF, a^SV^Vg , 

The tensors appearing above arise after the evaluation of the trace in (1) and the transformation of the fields to the 
basis with good quantum numbers. The connection between these coupling tensors and the original ones appearing in 
(1) which can be found in [27, 28] is given by: 

8 8 8 

i : j,k=0 i,j,k : l—0 i,j,k,l—0 

(3) 

The transformations preserve the symmetry structure of the tensors, that is G Q /3 7 and -F a /37<5 are completely symmetric 
and H a /3,-yS is symmetric upon the interchange of two indices which are on the same side of the comma. 

The tree-level mass square matrices are not diagonal in the 0,3,8 subspace, but since they are real and symmetric 
diagonalization is achieved with an orthogonal transformation. The tree-level orthogonal matrices in the scalar and 
pseudoscalar sectors are denoted with Os and Op, respectively. Denoting the eigenvalues of the pseudoscalar and the 
scalar 3x3 mass matrices in the 0,3,8 sector with Xp ,{ m i n ,mid,max} an d As^jn^j^max} the tree-level masses of the 
mesons are as follows 



m l+ = m l- = m P.12> m l+ = m l- = m S,12, 

m o = Ap m in, m 2 = As m id, 



K+ = m K- = TO P,45' m K + = m «- = m SA5> (4) 

2 2 2 2 2 2 

Tip 67, m K o — fn R a — m s,67i 



I)) 

m 2 = Ap, m id, ml = As jmin , 

m v , = Ap. max , m /o = ^S,mai 



Note, that some of the tree-level masses of scalars and pseudoscalars coincide. As we will see, the introduction of 
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TABLE I: The tj a coefficients appearing in the equations of state (6), (7), (8). The summation index 5 goes over 0,3,8. 



isospin and hypercharge chemical potentials distinguish between the particles and as a result all the one-loop pole 
masses will be different for /!/, [iy 0. 
The tree-level fermion masses are: 



9f 



(V2i> + V3v 3 + v 8 ), M d = -?—(V2vo - V3v 3 + v 8 ), 
12 v 12 



9f 



(V2v -2v s ). 



(5) 



The evolution of the condensates with the temperature or/and the chemical potentials is determined by the three 
equations of state 



a=l,2,4...7 
dcr 3 



7S{0,3,8} 



V3 



+ -^vo + (gi + y)( u 3 + v l) + (.9i + .92)^0 + V2g 2 v v^j v 3 - e 3 



E t3 fJflf a )-3 E [(O^SsOs) ii (a J a 1 ) + (OlP 3 P ) ii (n^)]+ 9 fN c ((uu)-(dd}), 



(7) 



/e{<r,7r} 
a=l,2,4...7 



7S{0,3,8} 



dL 



= (- — > = m 2 v 8 + -j=v a v 8 + -^=(w| - u l) + (.9i + tt)( u 8 + ^3)^8 + — h( v 3 ~ v l) v o + (ffi + 92)vqV 8 - e 8 



E tS f,MlU)-3 E [(Os^O s ) 77 ( ( T 7( T 7 ) + (0?P 8 Op) 77 (7r 7 7r 7 



(1=1,2,4.. .7 
.9F 



7£{0,3,8} 



2V3 



N c ((uu) + (dd) - 2(88)), 



(8) 



where in the mixing sector cr 7 stands for a, a[j,/o and similarly tt 7 denotes 7r°, 77, 77' as 7 = 0,3,8, respectively. 



denotes the antiparticle of / a , that is e.g. for / = a and a = 1 one has a\ = and cr{ 



In this notation 



(fa fa) = Tg(mf a ), (qq) = —4mqTp(m q ), where T^(mf a ) and Tp(m q ) stands for the bosonic, and the fermionic 
tadpole integrals, respectively. These integrals are given in Appendix B of [19]. The coefficients 1 1 a are listed in 
Table I. In the mixing sector, that is for 7 = 0,3, 8, the 3x3 matrices read: 



Sj — Go — -vqF^o — -v 3 F y3 — -vgFys, 

4 - 4 ~ 4 - 
P 7 = G + -v H 7 o + -v 3 H l3 + -v 8 H y8 , 



(9) 



with the definition: (G 7 ) Q /3 = G 7Q /j, {F 1 s) a p = F a p~ ( s, and (Hys)a/3 = H a p r/ g. All the indices run through 0, 3, or 8. 
A. One-loop parametrization of the model at zero temperature and density 



One has some freedom in choosing the set of equations which determines the 13 parameters of the model, namely 
m o 1 A j fi 1 9 1 9f , vo , v 3 , v 8 , e , £3 , e 8 and lf,lb- These latter two parameters are the fermionic and bosonic renormalization 
scales. For the parametrization we follow the method described in [19] where the renormalization of the model was 
also discussed. The only difference in the present case is the appearance of v 3 and £3. Since at zero temperature 
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FIG. 1: Diagrammatic representation of the one-loop pseudoscalar self-energies used for the parametrization. The label associated 
to the line denotes the propagating particle. 

and densities the effect of isospin breaking is small we use the same values for If and h as in [19] where these were 
determined by minimizing the deviation of the predicted mass spectrum from the physical one. The external fields are 
determined from the equations of state (6), (7), (8) once the remaining 8 parameters are known. 

In order to avoid the appearance of negative propagator mass squares in the one-loop finite temperature calculations 
in the broken symmetry phase we use the Optimized Perturbation Theory (OPT) of Ref. [29]. This amounts to replace 
the mass parameter — to 2 , in the Lagrangian with an effective, eventually temperature-dependent, mass parameter to 2 : 

L mass = ^m 2 TrM^M - -(ml + m 2 )Tr NlHl = -m 2 TrM^M - ^Am 2 TiM^M. (10) 

The counterterm Am 2 is taken into account first at one-loop level, while m 2 will replace to 2 , in all the tree-level masses 
and is determined using the criterion of fastest apparent convergence (FAC). We have chosen to implement this criterion 
by requiring that for ir + the one-loop mass calculated at vanishing external momentum stays equal to the tree-level mass 
(M n + = m n +). We have checked that imposing this equation for the neutral pion rather than the charged one results in 
no significant changes in the parameters. We note here that we were forced to use the definition M 2 + = —iG^ 1 (p = 0) 
instead of defining the one-loop mass as the pole of the propagator because in this latter case the solution to the gap 
equation, to be presented below, ceases to exist above a certain temperature, in accordance with previous investigations 
using the OPT [18, 29]. 

As described in details in [19], with the application of FAC one can eliminate the effective mass parameter to 2 in 
favor of the tree-level pion mass to 2 + in all the other tree-level masses of the propagators used to calculate the one-loop 
self-energies: 

™ 2 = ml + + -^U - -j=v a - ^g 2 v v s - (<?i + f ) v 2 Q - U + ^) v 2 3 ~ (fll + f ) < (11) 
In this way one obtains the following gap equation 

2 2 C C V^2 / gi \ 2 / 3(?2 \ n / <?2 \ 2 l 2 2/2\\ 

m«+ = ~m - -j=v + -j=v 8 + —g 2 voV8 + [gi + y J v + I 51 + — v 3+ [gi + — j + ReE 7r + [p = 0, m l (m n+ )), 

(12) 

where E ff + denotes the self-energy of ir + shown diagrammatically in Fig. 1. Equation (12) is the first from a set of four 
coupled non-linear equations which determines m 2 ,, fi, /2, g, if one knows gp, vq, V3, v&. Two further equations of the 
set are given by the one-loop equation for the 77 and K + pole masses 

M| = -to 2 



Op M( iee + ReE ,3, 8 (p 2 = Ml) P , (13) 
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r 2 2 c , C C 2 52 2^2 

- + = - m ° ~ vf' + i V3 vT 8 + 7E V0V3 ~ vT 0W8 + 71 



9i + f) vl + (, 91 + |) « 2 + (.91 + 7 -f) vl+ReX K+ (p 2 = M| + ), (14) 

where -M 2 ree is the tree-level mass squared matrix of the mixing sector without the mass parameter to 2 , the orthogonal 
matrix Op diagonalizes the expression in the round bracket, and So, 3, 8 is the self-energy matrix of the pseudoscalar 
mixing sector. This matrix is determined numerically. The last equation in the set is the FAC criterion for the kaon, 
which requires ReE^+(p 2 = M K+ ) — Am 2 = 0. 

The parameters gp,v ,v 3 ,v s are determined as follows. A linear combination of v and v s is determined by the 
tree-level PCAC relation for the pion decay constant (see Appendix of [28]) 

[2 1 

U:=d lla v a = ^-v + -^v 8 . (15) 

One can see from (5) that the same linear combination enters the expression of the average mass of the two light 
constituent quarks, so that the Yukawa coupling is given by gF = {M u + M d )/ f^. Another linear combination of vq 
and vs appears in the expression of M s in (5) which together with the PCAC relation (15) determines Vq and vs'- 

vo = n u { l + wrw d ) ' vs = w u i 1 - M^rk) ■ (16) 

The remaining parameter, V3 is obtained by requiring that the difference between the tree-level masses of tt + and 7r° 
equals the physical value (Am n ): 

-[OlMl ee O P ] n = (A mTr f. (17) 

This equation has two roots for V3, a negative and a positive one. The positive root would give m K a < m K + for the 
kaon masses. Since the opposite relation holds in nature, we choose the negative solution which is the physically valid 
one. 

We use the following values for the physical quantities: = 139.57 MeV, Am^ = 4.594 MeV, M K + = 493.677 MeV, 
Mr, = 547.8 MeV, f v = 93 MeV, (M u + M d )/2 = 313 MeV, M s = 530 MeV and in addition l b = 520 MeV and 
lj = 1210 MeV for the two renormalization scales. 

B. Introduction of the chemical potentials 

The introduction of the chemical potential for a system with a set of conserved charge operators is reviewed below. 
For vanishing external fields the Lagrangian (1) is invariant under the following global vector transformations 

M -> e - iaaG Me laaG = M - ia G [G, M] + 0{a 2 G ), 

tfj er lOLGG i) = V - ia G ip + 0(a 2 G ), ^ 

where G denotes the representation of the baryon (B), third component of the isospin (I) and hypercharge (Y) operators 
which are related to the diagonal generators as B = ^/|Ao, / = ^3 and Y = -^^8- The coefficients in front of the 
diagonal matrices are chosen such as to obtain the right quantum numbers when applying the operators on the quark 
fields. The consequence of this symmetry is the existence of conserved Noether vector-currents 



The conserved charge is defined as Q G = J d 3 xJ G (x). In terms of particle number operators the conserved baryon, 
isospin and hypercharges read as 

Q B = ^(Nu+Ni + N.-Ne-Ng-Ng), (20) 

Q 1 = \(N U - N, a -N d + N 3 + N K+ - N K - + N- K o - N K o + N K+ - N K - + N R0 - N K o) 

+N ai -K-+N, + -N^, (21) 

Q Y = ^ (N u - N a + N d - N d - 2N S + 2N S ) + N K+ - N K - + N K » - N R » + N K+ - N K - + N K o - N R0 . (22) 



c 



Note the different sign of Nj^a , A^o , N K o , N R o in Q 1 relative to Q Y . This is because the particles K° 



fall into 



different doublets from the point of view of the I3 and Y quantum numbers: K a ,K + and K~ ,K U form a .Z3 doublet 
while K°, K~ and K + ,K° form a Y doublet (likewise for scalars). 
The statistical density matrix of the system is given by 



exp[-/3(iJ-// G Q G )], 



(23) 



with G going over B,I,Y in the summation over this index. Using (20), (21), (22) one can rewrite (23) by regrouping 
the terms in the exponent according to different number operators and obtain p = exp[— [3(H — /LtjJVj)], where i goes 
over all the non-singlet particles to which the following chemical potentials were introduced in terms of /j,b, fJ-Y- 



Hu = - 


~Hu = 




1 


1 


= - 


-Md = 


1 

r s - 


1 

2 M/ + 


1 

3^, 


= - 


-Ps = 


1 


2 




M a + = 


= 


= 




- = fJ'i, 


= 

M«° = 


[i K Q = 


= 


1 

= -Ma-- = g 
1 

- -Ma° — — 2 



(24) 



The singlet particles (ir°, 77, rf , af], cr, and /o) do not contribute to the conserved charges and in consequence no 
chemical potential is introduced for them. By looking at (24) one can see that different members of a given multiplet 
(e.g. 7r + and 7r~) acquire a different combination of the baryon, isospin and hypercharge chemical potentials, which 
means that the chemical potentials remove completely the degeneracy between the members of the multiplets which we 
observe in the vacuum, both at tree and one-loop level. We have to keep track of the effect of 21 individually different 
particles, which makes things more complicated than in previous studies of this model. 

The effect of the chemical potentials is taken into account through the propagators which are introduced using 
the definition familiar from the theory of many-body systems. The relativistic formalism was developed in [30] and 
is reviewed in Appendix A, where the calculation of the self-energy using the finite-density Green's function is also 
sketched. 

In order to see explicitly that the particle and its antiparticle reflect differently the presence of a finite density medium 
we give here the tree- level propagators of K + and K ~ : 



G K+ (k) = 

G K - (k) - 
1 



2E k 

i 

2K 



l+n K +(E k ) n K +(E k ) l + n K -(E i! ) n K -(E k ) 



ko - Ek + ie 
1 + n K - (E k ) 
ko — Ek + ie 



ko — -Ek — ie ko + E^ — ie 
n K -(Ek) l + n K +(E k ) 



ko — -Ek 



ic 



ko + -E k + ie 

n K + (Ek) 
ko + -Ek + ie 



(25) 



and -Ep = 



p 2 + m 2 K± . The interpretation of the terms on the right hand side of 



where n K± (E p ) = - (E ^^ ±) _ 

(25) is as follows (from left to right): addition of a particle, removal of a particle, addition of an antiparticle, removal 
of an antiparticle. Note, that in the propagator of the K + the particle is K + and the antiparticle is K~ , while in the 
propagator of the K~ the particle is K~ and the antiparticle is K + . 

For all the other scalar and pseudoscalar fields the propagators can be written analogously using the chemical 
potentials defined in (24). For the fermions the propagators are given in Appendix A. 



III. THERMODYNAMICS OF THE MODEL AT FINITE DENSITY 



A. The influence of /i/ and /iy on the CEP 



With the parameters fixed in the previous section, we can solve the model at finite temperature and density using 
the formalism described in Section II B and in Appendix A. One calculates the 1-loop integrals entering the finite 
temperature and density version of the equations which determine the state of the system: the three equations of 
state (6), (7), (8) and the gap-equation for m 7r + (12). The relevant integrals are given in Appendix A. An observed 
smooth variation of the order parameters with the intensive parameter (T, or /j,b,i,y) indicates analytic crossover type 
transition. A first order phase transition is signaled by the multivaluedness of either one of the three condensates in a 
given range of variation of the intensive parameter. The point where by varying some parameter(s) the nature of the 



7 



> 

CD 



-0.05 
-0.055 

-0.06 
-0.065 

-0.07 
-0.075 

-0.08 
-0.085 

-0.09 
-0.095 
-0.1 





. 500 




- \. 800 
\ 900 ', 






950 




\ ref. curve 








987.9 \ / 






7 /*•■ 


'//*•-''' 


Hi=30 MeV ■ 









> 

CD 



50 



100 150 
T[MeV] 



200 



250 



-0.05 
-0.055 

-0.06 
-0.065 

-0.07 ; 
-0.075 : 

-0.08 
-0.085 

-0.09 







H B =500 MeV 




ref. curve 



.70/ 



50 



100 150 
T[MeV] 



200 



250 



FIG. 2: The generic temperature dependence of V3 for a crossover transition: on the l.h.s. /is changes and /ii = 70 MeV, on the 
r.h.s. fii changes while /is = 500 MeV. On both panels the reference curve (ref. curve) refers to the case (J,b = fJ>i = fJ>Y = 0. 

phase transition changes from crossover to first order one corresponds to a second order phase transition. 

The critical end point (CEP) is a second order phase transition point on the /is — T plane where by increasing fiB the 
phase transition as a function of T changes from crossover to first order (/ij and /zy are kept constant) . At vanishing /ij 
and fly the CEP is located in the point (T, hb)gep = (63.08, 960.8) MeV. The pseudo-critical temperature at vanishing 
chemical potentials is T c {^b,i,y = 0) = 157.98 MeV. 

Here it is important to note that, with the explicit isospin breaking taken into account, these values have significantly 
changed with respect to those obtained without isospin breaking at all (neither explicit nor spontaneous): (T, ^b)cep = 
(74.83,895.38) MeV and T c (hb,i,y = 0) = 154.84 MeV [19]. At first sight this is surprising since we have seen that at 
T = /j, = the effect of the explicit symmetry breaking is minimal. The difference is due to the behavior of the ^3 with 
the temperature. Without explicit isospin symmetry breaking 1*3 is identically zero for /ij = 0. When £3 7^ one can see 
by looking at the reference curve of Fig. 2 that with increasing temperature V3 is decreasing significantly compared to its 
T = value and reaches a minimum around the phase transition point where the influence of V3 becomes the strongest. 
The l.h.s. panel of Fig. 2 shows that the baryochemical potential magnifies this effect, implying that approaching the 
CEP the influence of i>3 is even stronger. According to our conjecture made in [19] that a smoother crossover at fiB = 
will require a larger value of /is to turn the phase-transition in T into a first order one, implying a larger value of 
Ms, cep, we can expect that the larger value of (1b, cep in the case of the explicit isospin breaking compared to the case 
in which the isospin breaking is absent corresponds to a higher value of the width of the chiral susceptibility AT c (xx)- 
Indeed, by looking at Fig. 3 one can see, that in the case with explicit isospin symmetry breaking AT c (xx) increased 
by ~ 20%, approaching the value of AT C (\^) = 28(5)(1) MeV at at [ij = 0. This value was obtained on the lattice 
in Ref. [31] after the extrapolation in the continuum limit was done, though in this lattice investigation the effect of 




30 40 
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FIG. 3: The width of the peak of the chiral susceptibility AT c (xx) as function of the isospin chemical potential with (without) 
explicit symmetry breaking external field £3 7^ (£3 = 0). In the chiral quark model \ = dx/de x where x = i/2/3(«o — va) is the 
non-strange condensate, e x = ^2/3(eo — £s) and as shown in [19] Xm ~ X X- 
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FIG. 4: The surfaces swept by the coordinates Tcep and 



Hb.cep [MeV] 
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b,cep of the critical end point as function of /ii and /iy- 



isospin breaking was not taken into account. It would be interesting to see whether similar effect is produced on lattice 
when m u ^ m&. 

Varying /i/ and [iy the location of the CEP in the [1b — T plane changes. Fig. 4 shows the surfaces swept by the 
two coordinates of the CEP as functions of /i/ and \iy . One can see that fXy has practically no influence on Tcep, 
which decreases very slowly, while with its increase «b,cep significantly decreases. The increase of [ii pushes the CEP 
towards higher values of /xb.cep and lower values of Tcep. This behavior is in concordance to what was previously 
written on the influence of v$ on the CEP at fXi = 0, since by looking at the left hand side of Fig. 2 one sees that at 
finite hi the isospin condensate v$ increases even more with the temperature. 

One can gain intuition on the way the chemical potentials fii an d influence the coordinates of the CEP by 
attempting a simple interpretation of our results in terms of generalized Clausius-Clapeyron equations applied to our 
system. The particle number and entropy densities of the two coexisting phases will be determined assuming an 
ideal gas of the quasiparticle degrees of freedom, which differ only in their respective masses on the two sides of the 
phase coexisting curves. The Clausius-Clapeyron equation successfully describe the slopes of phase coexistence curves 
of strong matter as functions of various chemical potentials and quark masses [11, 32, 33]. They are derived from 
the Gibbs-Duhem relation which connect the variation of the intensive thermodynamical parameters of a macroscopic 
system: 



dp = sdT + riBd[AB + nidfii + nydfiy. 



(26) 



Here ns,ny,n/ are the particle number densities and s is the entropy density. Keeping the pressure plus any other 
two intensive parameters constant one finds the following set of conditions for the phase coexistence when one varies 
the remaining two intensive parameters along the coexistence "surface": 



dT 



d\iB 



~As~' 



d[iB 



dT 



dfiY 



Any 



Any dfii 



Aub djii 



7>y 



dT 
dfj,i 



Ani 



An/ 



An B 



(27) 



On the right hand side of the equations above A refers to the difference of the values of a given extensive quantity 
in the symmetric and broken symmetry phase. In the two coexisting phases the relevant particle number and/or 
entropy densities (tig, G = B,I,Y, and s) can be calculated from the partition function using the formulas uq = 
TV~~ 1 dhi Z / due and s = V~ 1 d(T\nZ)/dT. Our simplified picture of the composition of the two phases in terms of 
non-interacting mixtures of 15 quasiparticles is given by 



kiZ = Vj2li(2si + 1)/ 7^ [fat + ln(l + a ie -^>-^) + ln(l + ai e~^ + ^) 



(28) 



where i £ tt , 7r°, K ± , K°, rj, rj' , u , u. , 
Si = for bosons, respectively. The energies u>. 



k° ,a, fo,u,d, s, 7i = N C7 at = 1, s; = 1/2 for fermions and 7, = oli = — 1, 
T.f are calculated with help of the tree-level mass expressions 



vV 



(4) after substituting into them the order parameter values determined in our field theoretical treatment for the two 
phases, that is by solving (6), (7), (8), and (12). 

The simple model predicts that Ang, Any and As is always positive, while An/ is always negative. Moreover, the 
following relations are obtained: Ang w Any, As > Ang and As > |An/|. The discontinuity of the particle number 
densities is determined by the contributions of essentially three quasiparticles: u 7 d 7 and n^ 1 . From our simple and 
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transparent model we get the sign and even the magnitude of the shifts of the CEP in agreement with Fig. 4 with the 
single exception of d^ij/dT . The ideal gas model does not reproduce the value of this derivative obtained by solving 
our model. We interpret this discrepancy as a result of the strong coupling between the (uu) ~ ^/2/3vo + \A/3 W 8 + V3 
and (dd) ~ y/2/3vo + y/T/3vg — U3 condensates not captured by the ideal gas approximation. As one can check also 
in [13] the strong coupling between these condensates reduce the temperature of the CEP when a finite fj,i is switched 
on. This is the same tendency we found in our field theoretical calculation. For the other three shifts it is the mass 
differences of the lightest quasiparticles of the effective model which exert the strongest influence on the variation of 
CEP position. 

B. Quasi- particle masses 

We turn to the study of the dependence of the tree-level masses and the one-loop pole masses on the temperature 
and the chemical potentials. The one-loop pole masses are determined as the zeros of the real part of the corresponding 
one-loop inverse propagators at vanishing spatial momentum. For example, the equation determining the one-loop 7r + 
mass reads: — ReG~+ (po = M v + , p = 0). If there are more than one solutions of this type of equations, then we 

follow that solution which in the vacuum lies closer to the physical mass. Usually this solution is lost as the temperature 
increases and some other solution is found. 




T [MeV] Hi [MeV] 

FIG. 5: The temperature dependence of the tree-level masses is shown in panels (a) and (b). The /i/ dependence of the one-loop 
pion masses for different values of fiB at T = (panel (c)) and T ^ (panel (d)). 

In Fig. 5(a) and Fig. 5(b) we see that the tree-level masses of ^ , tt° and a clearly reflect the restoration of the SU(2) 
symmetry at high temperature. This is not shown by the masses of Oq , a® and T). We cannot go to higher values of the 
temperature because at T ~ 252 MeV the non-strange condensate x decreases below the value of V3 and the tree-level 
mass of the u quark turns into negative. At this temperature there is still no sign for the tendency of the SU (3) chiral 
partners to become degenerate. 

In Fig. 5(c) we can see the dependence of the charged and neutral pion masses on the isospin chemical potential. 
The charged pions have by far the most significant dependence on /zj from all of the charged pseudoscalar mesons. At 
T = /j,b = the splitting between n + and ir~ is controlled by the bubble diagram involving tt + and n~ respectively 
(see Fig. 1) and the splitting point is at [ii ~ m n . One can see that at T = hb = the mass of tt° depends mildly on 
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111. This dependence intensifies with the increase of T and /is, but it remains true that the dependence on [ij is less 
strong than for the case of m 7T ± . It is interesting to note that for large values of fiB , when H u /d > m u/d and the fermion 
bubble contributes to the one-loop self-energies, the shape of the m n ±(fj,i) curves changes: m„+ starts to increase with 
Hi and the increase of m T - with \ij is slowed down and eventually turned over into a decrease in a given interval of Hi. 
Panel (d) shows that the increase of the temperature has a similar effect as /is in that it turns over the /^-dependence 
of ra n ± with respect to the behavior at T = \i& = starting at a low value of /xj. 
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FIG. 6: Estimation of the topological susceptibility through the Witten-Veneziano mass formula using one-loop and tree-level 
masses. The T-dependent pion decay constant (fn(T)) is approximated with the non-strange condensate x. The value of the 
pseudocritical temperature is T c = 157.98 MeV. In the inserted figure A//^ is plotted based on the tree and one-loop masses. 

In Fig. 6 we plot, both at tree and at one-loop level a combination of the masses and the pion decay constant 
A = (mjj + m 2 , — 2m|-)/^/6, which through the Witten-Veneziano mass formula [34, 35] 

^fx T = ™% + m%, - 2m 2 K (29) 

with Nf = 3, can be considered as an estimation of the topological susceptibility x T (T), which plays a crucial role in the 
phenomenology of the U(1)a anomaly (see e.g. [36, 37] for recent studies in terms of effective descriptions). In principle 
X T (T) can also be computed directly in our model if the quantity corresponding to the topological charge density Qt of 
the QCD is extracted. This can be done by comparing the four-divergence of the singlet axial vector current, which in 
QCD involves the U(1)a anomaly term with the corresponding current of the chiral quark model. Since the determinant 
term of Eq. (1) breaks the U(1)a symmetry, the correspondence is Qt ~ g(det(M) — det(M^)) = glmdet M. 

The decrease with T of the estimated Xt{T) seen in Fig. 6 doesn't mean the restoration of the U{1)a symmetry, since 
through fn{T), X T (T) is dominated by the restoration of the chiral symmetry. In view of (29) this can be also seen on 
the inserted figure of Fig. 6. However, the fact that at T = the estimated X T (T) is so close to the value obtained on 
the lattice in [38] and the curve itself stays within 10% of the lattice points, could imply that the effective restoration 
of the U(1)a symmetry, if contained in the lattice data 1 , could be implemented in an effective description based on the 
chiral quark model. In the NJL model the lattice result on Xt(^) [38] is converted into the temperature-dependence 
of the strength of the determinant term, by fitting it with the explicit formula of the susceptibility calculated in [39]. 

IV. CONCLUSIONS 

In this paper we studied the influence of the isospin and hypercharge chemical potentials on the [Ib — T chiral phase 
diagram of the three flavored chiral constituent quark model with explicitly broken SU(3)l x SU(3)r symmetry. The 
model was parametrized at one-loop level and optimized perturbation theory was used for the resummation of the 
perturbative series. Only one critical end point (CEP) is found for both spontaneous and explicit isospin breaking. In 
the latter case, based on the width of the peak of the chiral susceptibility, the crossover transition at fJ>B,i,Y = is 
found to be weaker than in the former case. Compared to the case without isospin breaking, in the case with explicit 
isospin breaking, the location of the CEP moves to a higher value of /is and a lower value of T. For fii = fj,y = 



1 Restoration of U(1)a symmetry requires that % T (T) decreases faster than f-^iT) with the increase of T so that X T (T)/f%(T) —> 0. 
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the coordinates of the CEP are: (T, hb)cep — (63.08,960.8) MeV. This value of ^b.cep is about three times larger 
than the value found on the lattice [40] and increases (decreases) linearly with fii (ply), while Tcep is two fifth of the 
lattice value and decreases slightly with the increase of \iy and significantly with the increase of Using an ideal 
gas picture and the generalized Clausius-Clapeyron equations we could interpret semiquantitatively with one exception 
the influence of [iy and fii chemical potentials on the CEP as resulting from the quasiparticle masses. We also studied 
the dependence of the charged and neutral one-loop pion masses on the isospin chemical potential at different values 
of the temperature and the baryon chemical potential. As a continuation of the present study, it would be interesting 
to investigate at what value of the fii do the charged pions condensate. 
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APPENDIX A: THE FORMALISM OF RELATIVISTIC MANY-BODY THEORY FOR A SYSTEM AT 

FINITE DENSITY AND TEMPERATURE 

We review below the method of relativistic many-body theory developed in [30] for the perturbative calculation of 
the self-energy at finite temperature and density. 

First we present the derivation of the tree-level Green's functions for K~ 7 K + , which depend both on the isospin 
and hypercharge chemical potentials. The field operators K~(x) and K + (x) are written in terms of creation and 
annihilation operators a + (p),6 + (p) and a(p),6(p), respectively, as 



A' + W = /|| 7 =(!. + (p) e "- + <.(P)^") 



Po=E p 

(Al) 

Pa=E p 



where E p = J p 2 + m 2 K± . This means that a + (p) creates a K + particle, b + (p) creates a K particle, etc. The 
operators have the usual non-zero commutators 

Hp), +(k)] = [&(p), b + (k)] = s(p - k). (A2) 

The two point functions for K~ and K + are defined as 



G K -(y - x) := {TK-{y)K+{x)) fi = Q(y Q - x ){K~ {y)K + {x)) p + Q(x - y )(K+(x)K-(y))p, 
G K+ {y~x) := (TK+(y)K-(x)) p = O(y - x Q )(K + (y)K-(x)) l3 + e(x - y )(K~(x)K + (y)) p , 

where the average is to be taken over a grand canonical ensemble, that is for an operator O one has 



(A3) 



(Oh - (A4) 

with H = H — \nQi. We make this distinction between K + and K~ propagators because the particle and its antiparticle 
feel differently the presence of the dense medium, resulting in a different mass dependence on the chemical potential. 
In our case this difference in the mass manifests itself first at one-loop level. 

Substituting (Al) into (A3), taking only the non-interacting part of the Hamiltonian H, with the help of the 
commutator relations given in (A2) and the Campbell-Baker-Hausdorff relation one evaluates the expectation values 
obtaining 

(a+(p)a(q)) p = S(p - q)n K+ (E p ), (6+(p)6(q)>/3 = S(p - q)n K+ (E p ), (A5) 

where n K +{E p ) = e!3{Ep _l K+ j _ 1 , n K +(E p ) = eP(Ep+ l K+ j _ 1 ■ Note, that n K +(E p ) = n K -(E p ). Using (A5) and 

the Fourier representation of 8(t) in (A3) one obtains in momentum space the K + and K~ propagators given in (25). 

Next, we calculate a one-loop bosonic bubble appearing in Fig. 1. With the standard rules of the perturbation theory, 
using the conventions of [41] the ir + self-energy is given by 

-iTi^+(y,x) = -4(3G 2 /3 7 + 4^2/3,7^5) (3Gi/3' 7 ' + ^Hip', J 'S'Vs')G. Kf}l . Kl3 {y,x)G rT _ iirT ^{y,x). (A6) 
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The first non-mixing bubble graph in the diagrammatic representation of given in Fig. 1 is obtained with the 
choice f3 = 4, = 5 implying 7 = 7, 7' = 6. Using that G156 = G247 and -ffi5,6<5 = H24.7S the contribution of this graph 
is 



K ~ (y, x) = y y~^ = -A 



K+ 



3G 24 7 + 4H 2 4,78VS 



G k +k- (y, x)G R o K o (y, x). 



(A7) 



The labels in the graph denote the field operators, e.g. on the left hand side ir creates a ir + particle. 
Going to momentum space one has 



(p) 



-4iV„ 



■K+K-K 



d^k 

—-^G K +(k)G R o(p- k) = 4V;+ K - K (.4(P, "»if + > 



(A8) 



where the vertex is V t +k-k° = 4 



V3 



and G K +(k) = G K + K -(k). 



Generally, at finite chemical potentials and temperature for a bosonic bubble diagram one calculates at vanishing 
spatial external momentum (p = 0) an integral of the form: 



l£(P0,TOi,/Ai,m 2 ,M2) 



d i k 



(27T) 

rf 3 k 1_ 

(2tt) 3 AE^E 2 



jG^Giip - k) 



p=0 



1 + m + n 2 



m - n 2 



m - n 2 



l + m + n 2 



Pa - Ei- E 2 p — Ei + E 2 p Q + E\ - E 2 p + E x + E 2 



,(A9) 



where for the propagators one uses a form similar to that in (25) and to arrive at the second equality one performs a 
contour integration in the complex energy plane. The distribution functions n, = (Ei) with Ei = v/k 5 " + mf contain 
the chemical potential for particle or antiparticle which is created by the fields of the vertex in the left hand side. 
We rewrite the integral (A9) as: 



/^(Po, mi,fxi,m 2 ,fj l 2) = Ib' T °(po,mi,m 2 ) 

1 2 



8ir 2 p 



dE\ E 2 - 



H(E) 



fh(E) 
poai - E poai + E 



(A10) 



where the remaining integral is evaluated numerically, V stands for principal value. The vacuum integral 
Ig T ~°(po,mi,m 2 ) is given by the expression (B4) of [19], i%i = l/(exp(/3(£ l — fa)) — 1) is the Bose-Einstein distribution 
and Oi = [l+ (-ly-^mj ~ m|)/jjg]/2. 

For fermions the method is identical to that used for the bosons. The fermion propagators for the constituent quarks 
u. u are defined as 



D u (y - x) := (Tu(y)u(x))p = 6(y - x )(u(y)u(x))p - 6(a; - y )(u(x)u(y))p, 
Du(y - x) := (Tu(y)u(x)} p = <d(y - x )(u(y)u(x))p - &{x - y ){u(x)u(y))p, 



(All) 



which in the momentum space read 



D u (k) = 
Du(k) = 



i(k ! + m u ) 
2El 

i( ft + m u ) 
2Eu 



W^k) 



H(Ev 



l-fu(Eu) 



fu(E u 



ka - E k + ie k — E k — ie k + E^ — ie k + E k + ie 



l-fu(E^) 



fu(E*) 



1 - /+(£k) 



H(Ev 



where f+(E p ) 



1 



ko - Ek + ie k — E k - ie k + _E k — ie fc + Ek + ie 
1 



(A12) 



and /„ (E p ) 



antiquarks. 

Then for the fermionic bubble appearing in the n + self-energy (see Fig. 1) one has 



are the distribution functions for u type quarks and 



d 4 fc 



9l 



l5 D d(k)D u (k +p) = ^N c I^(p,md,^ m u, /"«)• 



(A13) 



Similarly to equation (A9) in case of fermions we use the integral: 

I^(pQ,mi^i^ 2 ,m 2 ) = -iTr j ^Di(k)^D 2 (k + p) 

Jk 



p=0 
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+2(pf ) - {mi - m 2 Y) 



d 3 k 



d 3 k 

(2tt) 3 
1 



^rift + fi - 1) + ^rift + f-2 - 1) 

£j\ Ej 2 



(2tt) 3 ^E x E 2 



i - n - /a 4 



A" - / 2 + 



i - /r - n 



-2T 



fi,T=O t 



p - Ei- E-2 po - Ei + E 2 po+Et- E 2 po + E\ + E 2 



{mx)-2Tp (m 2 ) + 2(p - (mi - m 2 ) )I% (p , mi, m 2 ) 



+T^°(m 1 )+T^ v (m 2 ) 



pi - (mi - m 2 f 
4-K 2 p Q 



i=i 



dE\ E 2 



f+(E) f-(E) 



p Q a,i - E pocii + E_ 



(A14) 



where = l/(exp((3(E =p fn)) + 1) is the Fermi-Dirac distribution. 
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